Optimizing momentum space DMRG using quantum information entropy 
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In order to optimize the ordering of the lattice sites in the momentum space and quantum chem- 
istry versions of the density matrix renormalization group (DMRG) method we have studied the 
separability and entanglement of the target state for the 1-D Hubbard model and various molecules. 
By analyzing the behavior of von Neumann and Neumann-Renyi entropies we have found criteria 
that help to fasten convergence. A new initialization procedure has been developed which maxi- 
mizes the Kullback-Leibler entropy and extends the active space (AS) in a dynamical fashion. The 
dynamically extended active space (DEAS) procedure reduces significantly the effective system size 
during the first half sweep and accelerates the speed of convergence of momentum space DMRG 
and quantum chemistry DMRG to a great extent. The effect of lattice site ordering on the number 
of block states to be kept during the RG procedure is also investigated. 
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I. INTRODUCTION 

The density matrix renormalization group (DMRG) 
method 0, 01 has been widely used in recent years to 
study coupled fermionic and spin chain problems. Its ap- 
plication got a new momentum during the past few years 
when it was reformulated to models denned in momen- 
tum space 3 1 (MS-DMRG) or to quantum chemistry cal- 
culations [lp| (QC-DMRG). The properties of the Hub- 
bard model (3 and small diatomic molecules H U E3 
have been studied using these methods. However, several 
new technical problems are raised by these versions of the 
DMRG that have to be solved to increase the efficiency 
of the method and to stabilize its performance. 

A main difference between MS-DMRG and homo- 
geneous lattice models studied by standard real space 
DMRG with periodic boundary condition is that in the 
latter case each lattice site is equivalent and carries the 
same amount of information. In contrast to this, k- 
points or molecular orbitals lying closer to or farther 
away from the Fermi surface have different information 
content. The method is very sensitive to the ordering of 
the k-points or molecular orbitals and an optimal order- 
ing would have a major impact on the performance of 
MS-DMRG H S H 03 In fact the method can loose 
the target state and converge to a local minimum if an 
inappropriate ordering is used. It has also been found us- 
ing the dynamical block state selection (DBSS) approach 
0] that the same accuracy can be achieved with more or 
less block states, depending on the ordering. 

The density matrices of composite systems, the separa- 
bility of states and the nature of entan glem ent have been 
extensively studied [H El El El ElM El El E3 in 
the past few years. Since MS-DMRG represents a com- 
posite system with long-range interactions the results of 
quantum information theory can be used to understand 
the criteria of convergence of MS-DMRG. 

Another major feature which hindered the powerful 
application of MS-DMRG is the lack of the so called in- 
finite lattice method which in the real space version gen- 



erates a relatively accurate starting configuration for the 
so called finite lattice method. In the initialization pro- 
cedure proposed by Xiang 0] the environment blocks of 
various lengths were generated in advance of the finite 
lattice method. It is expected that a better procedure 
can be developed to generate the environment blocks by 
taking into account the change of the system block states 
during the process of renormalization. This is also crucial 
for reaching faster the crossover between the environment 
error and the truncation error. Then using the DBSS ap- 
proach the accuracy is controlled by the threshold value 
of the truncation error fixed in advance of the calcu- 
lations. The new procedure relies on finding the most 
important states carrying the largest information. The 
space of these states will be called active space (AS) fol- 
lowing a similar notation in quantum chemistry, namely 
the CAS of the complete active space self consistent field 
(CASSCF) method used to define the orbitals for the CI 
treatment. As it has been shown long ago |2fJ in the 
multi reference configuration interaction (MRCI) calcu- 
lations the convergence depends not only on the size of 
the active space but on other constraints of the numerical 
treatment as well. It is, therefore, a very important task 
to develop a protocol that extends the active space more 
effectively. 

From the point of view of synergetics DMRG can be 
interpreted as a dynamical system. In this analogy the 
response of a given model system to incident messages 
is studied and the change of the relative importance of 
messages is determined as the method converges to an 
attractor. Therefore, besides the practical importance to 
improve the MS-DMRG procedure, the study of the in- 
teraction of the subsystem blocks of DMRG as a function 
of the ordering of lattice sites is a very interesting ques- 
tion from the point of view of information theory, syner- 
getics and quantum data compression [2l], [22, [23, |24( ■ 

Our aim in this paper is 1) to study the criteria of con- 
vergence of MS-DMRG by analyzing the structure of the 
superblock and subsystem density matrices and quanti- 
ties used in quantum information theory, 2) to develop 
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a more efficient initialization procedure that collects the 
most important block states required to describe the to- 
tal system in a better way and 3) to study the effect of 
ordering on quantum data compression. 

The setup of the paper is as follows. In Sec. II we 
describe the theoretical background of separability and 
entanglement, mutual entropy and relative importance 
of incident messages, and generation of information en- 
tropy. Numerical investigation of these quantities in the 
context of MS-DMRG and QC-DMRG is presented in 
Sec. Ill and optimization of ordering is shown in Sec. IV. 
Sec. V. is devoted to the main steps of a new protocol 
that uses a dynamically extended active space to improve 
the initialization procedure. It is shown how the Abelian 
point-group symmetry can be used in the framework of 
QC-DMRG. Numerical results obtained for the 1-D Hub- 
bard model and for various molecules are presented in 
Sec. VI. and the effect of ordering on quantum compres- 
sion is also analyzed in some detail. The summary of 
our conclusions and future prospectives are presented in 
Sec. VII. 



II. THEORETICAL BACKGROUND 

A. Separability versus entanglement 

In general, if a finite system is divided into smaller 
subsystems (blocks) the Hamiltonian of the finite sys- 
tem is constructed from terms acting inside the blocks 
and the interaction terms among the blocks. The Hilbert 
space (G) of dimension N of the system is formed from 
the direct product states of the subsystem basis states. 
In particular, if Q describes a composite system with m 
subsystems then, 



G = ®?=xGi, dimQ = f[Ni = N. 



(1) 



In general, states of Q can be pure states or mixed states 
described by the density matrix written as 



(2) 



where \ipi) are eigenstates of the Hamiltonian acting on 
Q and YliPi = 1- The density matrix p has the following 
properties: (i) Trp = 1, (ii) p is a positive operator, i.e., 
Tr(pP) > for any projector P, (iii) p can be represented 
by its spectral decomposition as 



jv 



N 



P=^UJ n Pn, ^ U n = 1, U! n > 0, 



(3) 



n=l 



where P n form a complete set of orthogonal projectors. 
A state p is called separable if it can be written in the 
form 
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FIG. 1: Schematic plot of system and environment block of 
DMRG and related quantities described in the text. 



where p\ are states on Gi, thus the subsystems are either 
not correlated or their correlation is purely classical. 

In the density matrix renormalization group method 
proposed by S. R. White PJQ the total system called the 
superblock has two sites with qi and q r degrees of freedom 
between the left and right blocks Bi and B r of dimensions 
Mi and M r , respectively. Thus the superblock can be 
considered as a specially constructed composite system. 
Its configuration is shown in Fig. ^ 

The so called target state (|^t)) of the superblock sys- 
tem which is the eigenstate that one wants to calculate is 
formed from the direct product states of the blocks and 
the sites. It can be a coherent or incoherent superpo- 
sition of several eigenstates. In the first case one deals 
with a pure state, while the latter corresponds to a mixed 
state. In this paper we examine only the first case, i.e., 
\iPt) is chosen to be a pure state. 

If we combine the Bi* composite system to one sub- 
system and *B r to another one, then the so called su- 
perblock Hamiltonian for such a bi-partite system (in 
this paper bi-partite system refers to a composite system 
with tti — 2 and not the conventional term used in solid 
state physics to denote a crystal with specially arranged 
sublattices) consists of interaction terms determined in 
the blocks, denoted &sTIl and Hr and interaction terms 
between the blocks denoted as Hlr- The relative con- 
tribution of each term to the superblock energy can be 
measured as 



{Hi 



(^tIE.Wt) 



(5) 



where i = L,R, LR and X)i(^) = (J2i 

Since the target state is a pure state it follows from 
the Schmidt decomposition that for |^t) € G = Gl®Gr, 
with dim£ L = M L , dimGR = Mr, M l x Mr = N, 



r<min(Mx, ,Mr) 



(6) 



i=l 



where |e<) ® form a bi-orthogonal basis (ej|ej 



(4) (fi\f: 



- Sij , and < u>i < 1 with the condition 
^jW- = 1. If r > 1 and in the range of a state p there 
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exists a |\I/) such that 

A = (i-ip- 1 !*)- 1 > 



1 + maxi^j(uJiUjj) 



A = > maxw, 2 



(7) 
(8) 



then according to Ref. |15| p is inseparable. If p has an 
eigenvector |\&) corresponding to the eigenvalue A such 
that the conditions in Eqs. 0-© hold, then p is sep- 
arable. In our case for a pure state |\&) = \%t) an d 
A = 1. Using singular value decomposition to generate 
the Schmidt coefficients one can easily check if the condi- 
tions of Eqs. 0-© hold or not for a given target state. 

Even if the target state is a pure state the subsystems 
can be in a mixed state pl and pr. A quantitative char- 
acterization of the degree of mixtures is provided by the 
von Neumann entropy 



S{P) = -Tr(plnp) 
and the so called participation number 



R( P ) 



1 



(9) 



(10) 



The Neumann entropy is zero for a pure state and S 
!*] Y for a totally mixed states with pj — j^I where / 
denotes the identity matrix of dimension N. The par- 
ticipation number varies from unity for pure states to N 
to the totally mixed states and can be interpreted as an 
effective number of the states in the mixture. The Renyi 
entropy S q (p) = (ln(Trp 9 ))/(l — q) with q > 1 can also be 
used to measure how much a given state is mixed. The 
necessary conditions for a densit y m atrix to be separable 
have been worked out by Peres [25| and Horodecki et al 

m 

A fundamental concept related to inseparability and 
non-locality of quantum mechanics is the entanglement. 
A typical example of the maximally entangled pure state 
for two spin- 1/2 particles is 



= 7! m)± 



IT). 



(ii) 



In order to measure the degree of entanglement between 
the blocks of DMRG one can make use of the entan- 
glement monotone defined as E : p(Gl ® Qr) — > R+ 
with (i) E(p) = if p is separable, (ii) E is convex and 
(iii) E is non-increasing (on average) under LOCC (local 
quantum operations (quantum operations on the left or 
right block) or classical communications). A particular 
entanglement monotone, the entanglement of formation 
mUlEl is defined as 



E F {p) = min^p i S(Tr L |V'i)(V'i|)> 



(12) 



where the minimum is taken over all the possible realiza- 
tions of the state /?, TrL is a partial trace with respect 
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FIG. 2: The decomposition of the finite system to subsys- 
tems as a function of iteration steps corresponding to various 
partitionings. 



to the right block, S is the Neumann entropy defined in 
Eq. (JHJ. For a pure state Eq. (|T2"|> can be easily calculated. 
A mixed state is entangled if it cannot be represented as 
a mixture of factorizable pure states and there has been 
a great effort to determine a measure of entanglement for 
mixed states of a bipartite system [2!|,|3(],E1 ■ l n f ac t' ac ~ 
cording to Bennett [32l | all inseparable mixed states have 
non-zero " entanglement of formation" which means that 
a non-zero amount of pure entangled states arc needed 
to build them. 



B. Mutual entropy and Kullback-Leibler entropy 

In the DMRG procedure the complete Hilbert spaces 
of the blocks are truncated, thus they generate only a 
restricted subspace of the total Hilbert space. The vari- 
ous partitionings of the finite system containing L lattice 
sites to subsystem blocks are obtained by systematically 
changing the sizes of the left and right blocks (I) and (r), 
respectively, with I + 2 + r = L as shown in Fig. |2 In 
our implementation, described later in the text, the first 
iteration step corresponds to / = 1 and r = L — 2 — I and 
I is increased as long as I = L — 3. The renormalization 
procedure is used to obtain better configurations of the 
block states for a given partitioning yielding a more accu- 
rate description of the total system. In the forward sweep 
the left and right blocks are called the system and envi- 
ronment block, respectively. In the backward sweep the 
right block becomes the system block and the left block 
the environment block. In the DMRG method, the renor- 
malized states of the system block are selected from the 
eigenstates of the reduced density matrix of the system 
block having the largest eigenvalues. The reduced den- 
sity matrix is formed from the target state as ph = Tr^p 
and pn = Ti'l/3 for the left and right subsystems, respec- 
tively, and p = \iPt){4>t\ is determined by the diagonal- 
ization of the superblock Hamiltonian. It is thus evident 
that since the various possible partitionings of the sys- 
tem have a strong effect on the structure of Hl, 'Hr,'Hlr 
and on the density matrix of the superblock system, the 
structure of the reduced density matrix of the system 
block also depends on the environment block to a great 
extent. In this respect, states with largest eigenvalues of 
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the reduced density matrix of the system block can be 
considered as dominant states while states with smaller 
weights as recessive states. The term recessive is used 
to indicate that such a block state gives no considerable 
contribution to the superblock wave function with that 
given environment block, however, it is possible that it 
will provide a considerable contribution when it interacts 
with an environment block in a subsequent sweep of the 
finite lattice algorithm. The M m i n parameter introduced 
in Ref. ensures that these recessive states are also car- 
ried on during the sweeping procedure until they might 
become a dominant state. 

For a given system block various environment blocks 
can be constructed and analogously to the genetic algo- 
rithm we can treat them as different species of a popu- 
lation with different information content. Alternatively, 
we can think of the environment blocks as sources of mes- 
sages. A meaning can be attributed to a message if the 
response of the receiver, in our case the system block, 
is taken into account. Thus one can define the relative 
importance of message ojj as the eigenvalues of the re- 
duced density matrix. In other words, even if the en- 
vironment block contains all the states of the restricted 
Hilbert space defined on the r = L — I — 2 sites of the 
environment block, its information content can be very 
small if the system block is truncated so much that after 
taking into account the conservation of quantum num- 
bers the Hilbert space of the total system (superblock) 
is reduced drastically. The opposite treatment when the 
system block is considered as the source of messages and 
the environment block as a receiver works in the same 
way. Therefore, one needs a protocol to measure the 
conditional (mutual) information content of the blocks. 

Before proceeding further we recall a few definitions 
from classical information theory in order to describe 
Neumann entropy and its connection to quantum infor- 
mation theory. Let us assume that there are two events 
(not independent in general) described by two ensembles 
(X L , X R ) and two sets of probability distributions of the 
elements denoted as p(x L ) and p(x R ). We have used the 
labels L and R to indicate that the two events will later 
be related to the DMRG blocks. The Shannon entropy 
[33| for the two sets is defined as 

Mi 

S t = -5>0r})lnp(4), 5>(4) = 1 ( 13 ) 

J=l 3 

where i = L, R. It is worth to note that the entropy 
depends only on the probability distribution, p[x % j). The 
largest uncertainty of an event corresponds to the uni- 
form distribution with p(x l j) — 1/Mj where M, is the 
number of elements. If the state of the event is known 
exactly, then S — since p(Xj) = 1 and all p(xj'^j) — 
0. Following the notation of Haken [34| the expression 
Eq. (|13fl can be interpreted as an average of the quantity 
over fj 

3 



where fj = — Inpfej) and the weight is p(xj). In this 
respect fj is the information content of the state with 
index j and p(xj) is the probability or relative frequency. 
If under a different condition we find a new relative fre- 
quency p(x'j) then the change of information is 

Aj = lnp(x'j) ~ lnp( Xj ). (15) 

To obtain the mean change of information, we average 
Eq. i|15|) over the new distribution function p{x'-) and 
obtain the so called Kullback information gain as 

K(p(x'), P (x)) = "£p(x'j)Aj = £K4) ln ^3 ( 16 ) 

j j 3 ' 

where ^2jP(xj) = 1 and ^2jP(x'j) = 1. Eq. i|16|) has an 
important property, namely 

K(p(x'),p(x))>0. (17) 

If X L and X R are not independent of each other then 
the so called mutual information I(L, R) quantifies the 
correlation between the two events, i.e., it gives the in- 
formation about X L provided X R is known. I(L, R) is 
written as 

I(L,R) = S L -S(L\R) 

— Sl + Sr — Slr (18) 
= S R -S(R\L), 

where the conditional entropy S(L\R) is given by 

S(L\R) = Slr — Sr 

= -Ei^P(a£>ay)kp(ay,x£) (19) 

and the total entropy Slr is calculated from the joint 
probability distribution of the two events p(xj',x R ). It 
is clear from Eq. (|18l) that I(L, R) is symmetric under 
the interchange of X L and X R and zero if and only if 
the two events are completely uncorrelated, i.e., when 
p(x L ,x R ) =p(x L )p(x R ). 

The quantum analogue of the Shannon entropy is the 
Neumann entropy. In the context of DMRG it can be 
defined with p(x l j) = u>j where are the eigenvalues 
of the reduced density matrices of the subsystems. The 
two events can be related to the left and right subsys- 
tem blocks. Since the target state was chosen to be a 
pure state it corresponds to zero Neumann entropy with 
S LR = 0, thus I(L, R) = S L + S R . Recalling the Schmidt 
decomposition of Eq. © based on singular value decom- 
position one easily obtains that pl and pr have the same 
set of non-zero eigenvalues (u>i = u>f = u> R ) and 

Sl = Sr = ujj lawj. (20) 

3 

In what follows this will be denoted by S and called the 
mutual entropy. According to Eq. (|20|l , S is the measure 
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of entanglement formation. In order to describe the re- 
lationship between entanglement and mutual entropy we 
can use Wootters's interpretation j2^| obtained in quan- 
tum information theory. Within the context of DMRG it 
can be stated that for any pure target state the entangle- 
ment measures the amount of quantum information that 
must be exchanged between the DMRG blocks in order 
to create the target state. 

For a given system block we can generate various en- 
vironment blocks representing " different conditions" and 
the quantum analogue of Eq. called Kullback-Leibler 
entropy or quantum relative entropy is written as 

K(p\\a)) = Tr(plnp-plna), (21) 

where p and a are two reduced density matrices of the 
system blocks corresponding to two different environment 
blocks. 



C. Information generation and annihilation 

In the above treatment a so called static situation was 
considered where the DMRG system block has a given 
size with given lattice sites (orbitals) and we adjust the 
configuration space of the environment block. During the 
sweeping procedure of DMRG, however, the partitioning 
of lattice sites changes and the system behaves dynami- 
cally. In this case for a given incident message (given or- 
dering and configuration space of the environment block) 
the algorithm drives the system to an attractor defined 
by the target state. As it has long been known, in a 
dynamical system, different incident messages can give 
rise to the same attractor, in which case one can speak 
of redundancy of messages. On the other hand, it is 
possible that one incident message can lead to two dif- 
ferent attractors due to the fluctuation of the system or 
change of the intrinsic parameter of the system. This ef- 
fect doubles the original information. As an example, it 
was found earlier [9j that if the symmetry of the target 
state was not restricted, then depending on the DMRG 
parameters the same ordering of lattice sites and initial 
conditions sometimes gave rise to the S^ OT — compo- 
nent of the triplet state or to the singlet state. Clearly 
the relative importance of the message uij depends not 
only on the dynamical system but also on the tasks it 
must perform. In order to determine the values of ujj of 
incident messages we have to consider the links between 
a message and the attractor into which the dynamical 
system is driven after receipt of this message. Following 
again the notation of Haken a single message can drive 
the system via fluctuations into several different attrac- 
tors which may occur with branching ratios Mjk with 
Mjk = 1. Then the relative importance uij is defined 

as 

,. _Vr , .(i) _ \ - M 3k , .(i) r99 ^ 



where e — ► 0. If there are several systems coupled one 
after the other, then for instance in a two step procedure 
we obtain 

-i=E^ ) -i 1) = E4i ) 1 4 2 1 U: ) . (23) 

fc fcl,fc2 

It is worth to emphasize that the recursion from u>( n > 
to to may depend on the path, namely on the ordering of 
lattice sites or molecule orbitals. In such a way we obtain 
an interference of messages and the relative importance of 
messages depends on the messages previously delivered to 
the receiver. In the general case, the relative importance 
of a messages will depend in a non-commuting way on the 
sequence of the messages. In this way the receiver (in our 
case the system block) is transformed by messages again 
and again and clearly the relative importance of messages 
will be a function of the iteration step. 

The meaning of Eq. (12;it is that with a given task or 
ensemble of tasks, this algorithm allows us to select the 
message to be sent, namely the one with the biggest uij. 
Using the relative importance of messages one can inves- 
tigate the information entropy of subsequent states of the 
system 

S<°) = -^c^ln^-, =-^w fc lnw fe , (24) 

3 fc 

whether a dynamical system annihilates, conserves or 
generates information. If ujk = 1 and Y] j luj = 1 

there is annihilation of information if S*' 1 ' < S^ -* or gen- 
eration of information if > S^K A dynamical system 
can be "sensitive" or "insensitive" for a given message. 
As it was described by Haken 0| it is an interesting 
problem to determine the minimum number of bits re- 
quired to realize a given attractor or to realize a give 
value of "relative importance". The M m i n parameter in- 
troduced in a previous work |!j had the same meaning. 
It is, however, a major task to determine those incident 
messages which have the largest uij in advance of the 
calculations. 



III. NUMERICAL INVESTIGATION OF 
NEUMANN ENTROPY AND ENTANGLEMENT 

The elements of information theory outlined in the pre- 
vious sections have been used in the DMRG studies of 
systems described by the Hamiltonian, 

n = J2 T V c l* c j<?+ E V ijkl4a C ]o-' C kT'Cla, (25) 
ija ijklaa' 

where denotes the matrix elements of the one-particle 
Hamiltonian and Vijki stands for the matrix elements 
of the electron interaction operator. Depending on the 
structure of Ty and Vijki this Hamiltonian can describe a 
molecule or a usual fermionic model in solid state physics, 
e.g., the Hubbard or extended Hubbard models in one or 
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Orbitals ordered accoding to the legend 

FIG. 3: Decay of orbital occupation number obtained in the 
full CI limit and the corresponding ordering of orbitals for the 
molecules studied in the paper. 



higher dimensions or for example coupled fermion chains. 
In the former case a one-dimensional chain is built up 
from the molecular orbitals that were obtained, e.g., in a 
suitable mean-field or MCSCF calculation and in the rest 
of the paper we use the numbering of orbitals correspond- 
ing to the output of these calculations. In our solid state 
physics applications the indices k, I denote momenta 
with ki = (2im)/L,—L/2 < n < L/2. For the one- 
dimensional Hubbard model 7y = —2tcos(ki)8(i — j) 
and Vijki = (U/L)S(i + j — k — I). In what follows 
U is given in units of t with t = 1. In the rest of 
the paper Hartree-Fock (HF) orbitals denote filled ki 
points between the Fermi surfaces (±fcp) in the limit of 
Vijki = where we use "Fermi surface" to denote sites 
where the occupation number of the sites drops to zeros 
for Vijki = 0. The full CI (FCI) energy is the exact solu- 
tion of Eq. H25|) for a finite chain with length L and for 
a given number of electrons with up and down spins. 



A. Analysis on small molecules 

First we have analyzed the density matrices and the 
energy of the subsystem blocks and the interaction be- 
tween the subsystem blocks by carrying out test calcula- 
tions on a very small system, namely the CH2 molecule 
by correlating 6 electrons on 13 orbitals. The system is 
so small, that it is hardly expected that DMRG would 
not converge to the attractor determined by the target 
state. In the ordering we relied on the occupation num- 
ber of the orbitals obtained in the full CI calculation by 
the MOLPRO program package j3f|. This is plotted in 
Fig.[3]for a few selected test molecules used in the present 
paper. The legend shows the corresponding orbitals with 
the original indices. It is worth to mention that for the 
CH2 molecule orbitals 1 and 10 are almost doubly occu- 
pied while 2 and 8 are almost occupied with up or down 
spins since the ground state is a triplet state. 

In a first attempt to generate a good environ- 
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Iteration step 

FIG. 4: Example for an ordering when relevant information 
of the system is collected into one DMRG block giving rise to 
vanishing interblock interactions and very low entropy indi- 
cating high separability of the target state and lack of quan- 
tum information exchange between the blocks. 



ment block, we have put the orbitals with the 
largest occupation numbers to the right end of the 
chain. This corresponds to the following ordering: 
[7, 12, 4, 9, 6, 5, 11, 13, 3, 8, 2, 10, 1]. Our result obtained 
for a truncation error TRE max = 10~ 8 is shown in Fig. 0] 
In subsequent panels we show the number of selected 
block states, the relative error of the energy, the rela- 
tive energy of the blocks and the interaction term, the 
mutual entropy, the site and block entropies and the 
site participation numbers. It is evident from the fig- 
ure, that although depending on TRE max the number of 
selected block states fluctuated, the energy did not con- 
verge to the FCI value. It was trapped at some local 
minimum. This extreme situation has often been found 
in our previous studies. Investigating the third panel of 
the figure, it is obvious that the right block alone pro- 
vides all the contribution to the superblock Hamiltonian 
as ((Hr) ~ — 1) and the interaction between the blocks 
vanishes ((Hlr) — 0), except at the turning points. Due 
to the lack of interaction of the blocks the mutual infor- 
mation entropy of the blocks remained close to zero thus 
no quantum information exchange was generated dur- 
ing the sweeping procedure as can be seen on the fourth 
panel of the figure. In order to analyze the information 
content of the subsystems we have also plotted in the fifth 
panel the calculated entropies of the left and right blocks 
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(Si and S r ) and the entropies of the two intermediate 
sites (S q i and S qr ). The sixth panel shows the partic- 
ipation numbers of the two intermediate sites (R q i and 
R qr ). It is clearly seen that sites lying close to the " Fermi 
surface" have much larger Neumann entropy (S q i,S qr ) 
and larger participation number (R q i,R qr ) than orbitals 
where the occupation number is close to zero or two. 
(This is determined by identifying the data points cor- 
responding to the ordering and the partitioning of the 
composite system for that given iteration step as was 
shown in Fig- El ) The meaning of site entropy and par- 
ticipation number can be explained very easily. In our 
case the maximally mixed state of an orbital described 
by the qi = q r = 4 (0, |, !, t j) basis states would corre- 
spond to S = In 4 and R — 4. For those sites which lie 
energetically far above ( below) the Fermi surface the |0) 
(| tl)) basis state in the reduced density matrix of the 
site would appear with very large probability and the 
remaining three basis states with vanishing probability 
giving rise to S ~ and R ~ 1. In contrast to this, 
the reduced density matrices of sites lying close to the 
Fermi surface will have a more uniform eigenvalue spec- 
trum corresponding to a finite value of the entropy and 
a participation number larger than unity. Analyzing the 
situation close to the turning points, we have to recall 
that in this case the 'Br subsystem contains two lattice 
sites, but orbitals 1,10,2,8 are very important compo- 
nents of the wave function. Therefore, both subsystem 
contain important states and a finite quantum informa- 
tion exchange between the blocks is necessary in order 
to generate the target state. As a test we have flipped 
the ordering and found a similar behavior, except that 
the peaks of the interaction terms and the entropies were 
shifted to the other end of the chain. 

As a next step, the interaction between the 
blocks was maximized by putting the orbitals al- 
ternatively to the two ends with decreasing occu- 
pation number. This gives the following ordering: 
[1,2,3,11,6,4,7,12,9,5,13,8,10]. Our result is shown 
in Fig. [5J It can be seen that the convergence became 
very fast and within one and a half sweeps (26 itera- 
tions steps) the error margin set by TRE max was reached. 
This means that the environment error was reduced sig- 
nificantly within one full sweep. Investigating the third 
panel one can see that there is always a strong interac- 
tion between the blocks as (Hlr) — 0.5 and both blocks 
provide equal amount of energy for the total system as 
(TLl) — (Hr) — —0.7. Analyzing the separability and 
entanglement of the target state the fourth panel shows 
that the system possesses large mutual entropy indicat- 
ing a large entanglement and a large amount of quantum 
information exchange between the blocks. Another ma- 
jor difference compared to the results obtained for the 
previous ordering is that the block entropies (Si and S r ) 
are always very large as can be seen in the fifth panel. 

Investigating the first panel of Fig. [S] one can see 
that the maximum number of block states selected dy- 
namically were often close to 250. Based on previ- 
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FIG. 5: Same as Fig. [I] but for an ordering when relevant 
information of the system is shared between DMRG blocks 
giving rise to strong interblock interaction and large quantum 
information exchange between the blocks. 



ous studies we expect that the Cuthill-McKee algo- 
rithm |36| gives a better ordering in the sense that a 
smaller number of block states have to be selected to 
achieve the same accuracy. Our result for the ordering 
[7, 6, 5, 4, 3, 1, 2, 9, 8, 13, 12, 10, 11] is shown in Fig. @ It 
can be seen that in fact the same accuracy was achieved 
as before with a smaller subset of block states (M max < 
130). At the same time the mutual entropy had usually 
large peaks when the energy improved significantly. An- 
alyzing the fifth and sixth panels one sees that orbitals 
which have large entropy are now at the center of the 
chain and the site entropies had maxima where the mu- 
tual entropy had sharp peaks. In addition, compared to 
the two previous orderings it is evident that the mini- 
mum of R was at larger values. Looking at the third 
and fourth panel we can find an extended region with 
large interblock couplings but with small entropy. For 
this region our result seems to indicate smaller entan- 
glement between the blocks thus it might correspond to 
a more pronounced finite classical interblock interaction 
with less quantum correlation between the blocks. Again 
the significant improvement in the energy happens when 
the two blocks have the same amount of intrablock en- 
ergy with a large quantum information exchange between 
the blocks. 
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FIG. 6: Same as Fig. Hand Fig.EJbut for the Cuthill-McKee 
ordering. 
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FIG. 7: Result for F2 molecule correlating 18 electrons on 18 
orbitals using Cuthill-McKee ordering. 



B. Analysis on larger molecules 

In order to search for an optimal ordering method we 
studied the quantities defined above for the molecules 
CH2, H2O, F2 and N2. We have found that in general the 
Cuthill-McKee ordering is not the best one. In the case 
of F2 molecule with 18 electrons on 18 orbitals DMRG 
did not converge at all, as can be seen in Fig. It is 
evident from the figure that the orbitals with larger par- 
ticipation number and larger entropy appear close to the 
two ends of the chain and the method performed rather 
badly. A large number of block states were selected (we 
have cut M max at 2000) and the result did not converge. 
Cuthill-McKee ordering resulted in a very similar dis- 
tribution of site entropies as in the second example for 
the CH 2 molecule. We have tested this kind of entropy 
distribution for the other molecules as well using larger 
basis sets and obtained very slow convergence or lack of 
convergence. Therefore, the good convergence obtained 
in the second example for CH2 molecule was due to the 
short chain length. 



C. Analysis on the 1-D Hubbard chain 

We have done similar calculations on the 1-D Hubbard 
model for various band fillings, chain lengths and U val- 
ues. Fig. [8] shows six different orderings of the k values 
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FIG. 8: Example for various orderings of ki points for the 
half-filled Hubbard chain with L=14. Doubly filled sites in 
the HF limit in the natural ordering are 4,5,6,7,8,9,10. 



for a half-filled L = 14 chain. Using the natural ordering, 
where the values — tt < ki < it follow each other with in- 
creasing ki, the doubly filled states of the HF state are on 
the sites 4,5,6,7,8,9,10, the others are empty. The states 
at the Fermi energy are on the sites 4,10. 

The first panel corresponds to the original ordering, 
while in the second panel the ordering of sites with neg- 
ative and positive momenta were reversed. In the third 
panel sites counted from the Fermi surface were ordered 
to one side of the chain similarly to the first example 
used for CH 2 . In the fourth ordering the Fermi points 
were moved to the two ends of the chain similarly to the 
second example used for CH 2 . The fifth and sixth pan- 
els show our attempts to place the sites near the Fermi 
surface to the center of the chain as in the third example 
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FIG. 9: Site entropy (S q ) obtained for the half-filled Hubbard 
chain with L = 14 using ordering shown on the 6th panel with 
U=0.5,l,2,4. 



for the CH2 molecule. 

We have found similar results as for the molecules, 
namely the third ordering gave the worst result. For 
small U values (U ~ 0.5, 1) the method did not converge 
even if M m i n = 800 was used. The fourth ordering gave 
very large mutual entropy and although a very large sub- 
set of block states were selected, the convergence was still 
very slow. It required six to seven sweeps for a half-filled 
chain with L = 14, U = 1, and TRE max = 10~ 6 . The 
best performance was obtained for the fifth and sixth or- 
derings. The first two orderings gave stable but consid- 
erably slower convergence. We have carried out several 
calculations with chain lengths up to 30 sites for various 
fillings and U values and found again that good conver- 
gence is obtained if the block entropies are larger then 
a critical value depending on the various model parame- 
ters. As an example the distribution of site entropy for 
the half-filled Hubbard model with L = 14 using the sixth 
ordering is shown in Fig. [5] As can be seen larger U val- 
ues gave larger values of site entropies and participation 
numbers which means more mixed states. We have also 
found that the best ordering is a function of the filling. 
This is because the algorithm will converge only if there 
is a finite interaction between the blocks for several steps, 
and therefore both blocks must contain sites with large 
entropies. At the same time the reduction of block en- 
tropies and the number of block states can be achieved 
by moving sites with largest entropies close to each other 
and closer to the center of the chain. As an example, 
for larger U = 2,4,20 values, where sites lying farther 
away from the Fermi surface become more important, 



we have found that DMRG blocks had finite interaction 
strength for all iteration steps but of course the number 
of maximally selected block states has increased signifi- 
cantly. As an indication for an L = 14 half-filled chain a 
10~ 5 absolute error of the energy has been reached with 
M max ~ 400 — 500 for U = 0.5, while the system selected 
out dynamically more then 2000 — 2500 states to reach an 
accuracy of 10 -4 for U — 4. For L — 18 these numbers 
are 600 - 700 and 3000 - 3500, respectively. This fact 
questions the efficiency of MS-DMRG when U is compa- 
rable to or larger than the band width. We expect that 
for the large U limit real space DMRG should provide 
more accurate results. 



IV. OPTIMIZING ORDERING 

It is very important to emphasize that according to 
Eq. (I23|) the dynamics of the system depends very much 
on the path it evolves, thus on the way it receives incident 
messages and the renormalized block states are formed. 
In DMRG the initial wave function does not point in the 
direction of the attractor and it is rotated systematically 
in the multidimensional space by the transformation ma- 
trices during each renormalization step. The sequence of 
rotation, however, is not commuting, therefore, as it has 
been seen in the calculations presented above the order- 
ing has a great impact on the structure of renormalized 
states. For fermionic models in solid state physics the op- 
timal ordering following from our earlier considerations 
is the one shown on the sixth panel in Fig. [S] It is a 
far more complicated task to find it in quantum chem- 
istry. In order to search for an optimal ordering we have 
carried out a kind of brute force calculation for the F2 
molecule. We have run QC-DMRG using 64 block states 
for one half sweep, then permuted two orbitals and moni- 
tored the obtained energy values. If an ordering produced 
a better result it was kept and a next ordering was ob- 
tained by permuting again two sites. For the F2 molecule 
one calculation took some eight seconds CPU time and 
more than 50 000 permutations were carried out. Our 
result for the best ordering is shown in Fig.^J It can be 
seen that the sites with largest entropy are close to the 
center of the chain but some sites with finite entropy are 
close to the two ends of the chain. This configuration is 
in agreement with what has been found for the Hubbard 
chain. It provides finite interaction between the blocks 
for several steps but the largest information exchange be- 
tween the blocks occurs when the superblock composite 
system is close to a symmetric configuration. We have 
investigated the other test molecules as well and obtained 
similar results. It is worth to mention, that for a given 
molecule the ordering which produced the best perfor- 
mance at the end of the first half sweep also gave the 
fastest convergence using more sweeps. 

As a summary of our result, we have obtained the fol- 
lowing rule to generate a good ordering, which of course 
could still be not the optimal one: 
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the target state became non-separable. Since the subsys- 
tem entropies and related quantities depend on the other 
subsystems, this suggests that a better initialization pro- 
cedure can be obtained if one starts with a larger mixture 
by increasing the block entropies. This also leads to the 
decrease of the volume of separable basis states of the 
target state. We have also found that the entropy of the 
blocks must be larger than a critical value. 

The standard initialization procedure outlined by Xi- 
ang, when subsystem blocks of various lengths (describ- 
ing various partitions) are generated in advance of the 
finite lattice algorithm does not guarantee this in general 
for the first sweep. Therefore, it is crucial to develop 
a new method which constructs the environment blocks 
by taking into account the change of the renormalized 
system block basis states, and which increases the block 
entropies above a threshold value during the initializa- 
tion procedure. In this section we present a procedure to 
generate the system and environment blocks giving large 
LOj eigenvalues, i.e., block entropies. 



FIG. 10: Result for F2 molecule correlating 18 electrons on 
18 orbitals using our numerically optimized ordering. 



A. The DEAS procedure 



Step 1. The Cuthill-McKee ordering is used and a 
DMRG calculation with some 100 block states is run and 
the entropies, participation numbers and orbital occupa- 
tion numbers are determined. 

Step 2. The orbitals with largest entropies (largest 
participation numbers, closest to the Fermi surface ) are 
moved to the center of the chain and the remaining or- 
bitals are distributed along the chain with decreasing en- 
tropies counted from the two ends of the chain. 

This method guarantees a finite interblock coupling for 
several steps but also reduces the number of block states. 
It is worth to mention, that if q r < M™ ln or q< < M™ ln 
then the sweeping procedure can be turned back, since 
the remaining part of the sweeping will not improve the 
environment. All calculations in the rest of the paper 
were obtained following this rule. 

V. THE INITIALIZATION PROCEDURE 

Besides ordering, the optimal performance of DMRG 
is strongly effected by the initial conditions or in other 
words by the initial block configurations. Looking at 
Fig. one can see that in the first three iteration steps 
the mutual entropy S was exactly zero, R q i and R qr was 
unity (i? — 1 = 0) while in the second sweep at the same 
partitioning their values became larger indicating a bet- 
ter environment block. In addition, for the first three 
iteration steps the eigenvalue spectrum of the reduced 
density matrix of the system block (pl) had one eigen- 
value of unity which according to Eq. JSJl means that 
the superblock eigenstate was separable. For the second 
sweep all eigenvalues were less than unity indicating that 



Let us consider first the partitioning where the left 
block (system block) contains one lattice site and the 
right block (environment) L — 3 lattice sites as was shown 
in Fig[3 The system block contains Mi = q states and 
the right block contains M r states when the conservation 
laws of total quantum numbers are taken into account. 
The various configurations of the right block interact dif- 
ferently with the left block. Since the density matrix 
of the system block depends strongly on the interaction 
between the two blocks the initial configurations of the 
right block have a strong effect on the renormalized block 
states of the left block. In order to guarantee the conver- 
gence of DMRG the mutual entropy should possess a fi- 
nite value, thus one needs a method which maximizes the 
Kullback-Leibler entropy given by Eq. (|21|) . i.e., a proto- 
col to optimize the right block configuration in advance 
of the calculation to have largest Neumann entropy. 

It is known for the Hubbard model, that for small U all 
the interesting physical properties are determined by sites 
that lie close to the ±kp Fermi points. They give also 
the largest contribution to the correlation energy. The 
most important excitations are shown in Fig. 1111 As the 
strength of U is increased this region gets larger and sites 
further away from the Fermi points will also contribute 
to the correlation energy and other physical properties. 
We can thus expect that the basis states formed from 
the sites around ±kp will correspond to a larger mixture 
(several ujj will be large) and we can form the incident 
message from these states. In order to achieve this we 
define a so called AS vector, a concept inherited from 
quantum chemistry which simply contains the most im- 
portant sites in a descending order with respect to their 
expected importance. For the Hubbard chain this vector 
is AS = [fci?,— kp,kp + l,— kp — l,kp — 1,— fcp + 1, ...]. In 
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FIG. 11: Dispersion relation of the 1-D Hubbard chain. Ar- 
rows indicate the most probable electron-hole excitations with 
small momentum transfer in the small U limit. 

QC-DMRG the AS vector is constructed in a self consis- 
tent way. First the AS vector is defined to include only 
the HF orbitals and a quick calculation is performed with 
some 100 block states using one or two sweeps. Next the 
AS vector is defined by taking orbitals with largest values 
of site entropies in a descending order. 

We have confirmed numerically that the entropy of the 
environment block (S r ) can be increased by using sites 
lying around the Fermi points. It is, however, not suffi- 
cient to increase S r only because it could happen that a 
block configuration with large S r could lead to vanishing 
mutual entropy (S) due to the lack of interaction be- 
tween the blocks. Once a good AS vector is established, 
the question arises how to construct the configuration 
space, i.e., how to fill the sites given by the elements 
of the AS vector. Therefore, one needs a protocol that 
extends optimally the active space in order to achieve 
fast convergence. In the following we describe the de- 
tails of the dynamically extended active space (DEAS) 
procedure which takes into account the change of the 
partitioning of the DMRG superblock and the renormal- 
ization procedure of the system block. 

Let us assume that we have a given ordering and we 
have a starting partitioning with I = 1 and r — L — 2 — 1 
as is shown in Fig. [5] 

Step 1. An AS vector is defined and another vector 
(HF vector) is defined which contains the Hartree-Fock 
orbitals. 

Step 2. Based on the ordering the AS and HF vectors 
are reordered and for the given partitioning the elements 
of the AS and HF vectors belonging to the right block 
are determined. 

Step 3. The q 3 bases states of the Bi • • subsystem 
are formed and their quantum numbers are determined. 

Step 4. On the site given by the first element of the 
AS vector the basis states (0, |, 4, T4) are formed and the 
HF sites are filled according to the HF configuration ex- 
cept for this site if it is in the HF vector. An example 
is shown in Fig. ^| for the half- filled Hubbard chain us- 
ing the original ordering and M™ m = 8. The quantum 
numbers of the obtained states are determined and those 
which have a matching component with any of the Bi • • 
basis states satisfying the conservation of total quantum 
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FIG. 12: Example for a right block configuration space ob- 
tained by the DEAS procedure for the half-filled Hubbard 
chain with L = 14, iV T = ATj. = 7, M r min = 8, I = 1, 
r = 11 and HF vector = [4,5,6,7,8,9,10] and AS vector 
= [4,10,3,11,5,9]. 



numbers are kept. Next the first two elements of the 
AS vector are taken and the same procedure is repeated 
until M™ m many states are selected. In this way one 
makes sure that states of the environment block have a 
matching component with the system block, thus they 
all produce contribution to the mutual entropy. 

Step 5. One iteration step of the standard DMRG 
procedure is carried out. In the next step then procedures 
of steps 2-5 are repeated by forming again the MiXqXq 
states of the B\ subsystem and the states of the right 
block matching the left block are determined. 

Step 6. Close to the turning points it can happen that 
for the selected states M r < M™ m thus the excluded 
states could never be recovered during later steps of the 
method. In order to avoid such problem all q r states are 
selected for the environment block if M r < M™ m . 

This procedure ensures that the size of the active space 
is extended dynamically as the system evolves and the 
left block is correlated from the very beginning with the 
most important states defined by the AS vector. In ad- 
dition, as the size of the left and right blocks changes 
the second step of the procedure adjusts the AS vector 
so that only the those elements are taken into account 
which belong to that given right block. 

As an example, Fig. ll3l shows the density matrix spec- 
tra, the calculated mutual entropy (S), environment 
block entropy (S r ) and Kullback-Leibler entropy (K) for 
three specially constructed right blocks with M™" = 64 
for the half- filled Hubbard chain with U = 1, L = 14 
using the original ordering. In the first example we have 
formed the right block basis states from the sites ly- 
ing close to the right hand side of the chain like in the 
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FIG. 13: The reduced density matrix eigenvalue spectrum, 
the mutual Neumann entropy (S), environment block entropy 
(S r ) and Kullback-Leibler (K) entropy calculated for three 
specially constructed environment blocks for the half-filled 
Hubbard chain with U = 1 and L = 14 using the original 
ordering. 



standard initialization procedure. The second example 
corresponds to a similar construction but the Hartree- 
Fock sites were doubly filled. In the third example the 
right block was set up according to the method described 
above. It is evident from the figure, that in the first ex- 
ample one eigenvalue of the reduced density matrix was 
one giving rise to zero mutual and environment block en- 
tropies and separability of the superblock eigen state. 
The second example corresponds to a finite but very 
small entropy while the entropy is increased for the third 
example. It is thus clear that a better environment block 
can be constructed by using the AS vector. 



B. Reducing the effective chain length during the 
initialization procedure 

Once the M r states of the environment block are found, 
taking into account the conservation laws of the quantum 
numbers, many sites will remain unfilled as can be seen 
on the example shown in Fig. If the right block con- 
figuration space is represented by an M r x r matrix then 
those columns (sites) which contain only the empty state 
(labeled as E) are not active, since the action of creation 
and annihilation operators on such empty sites gives no 
matrix elements. Therefore the i,j,k,l summation in- 
dices of Eq. (|25[l can be restricted only to the remaining 
non-empty sites (NE) giving rise to a much smaller effec- 
tive size of the chain. 

The non-empty sites fall into two further categories. 
Those columns (sites) which contain only doubly filled 
sites (D) and those called active sites (A) which con- 
tain various basis states (0, |, f , If). Several steps of the 
DMRG matrix algebra can be restricted to the active 



sites only and various restrictions are also given if an in- 
dex corresponds to a doubly filled site. For example the 
matrix algebra for quadratic auxiliary operators having 
the form CiaCjp with a, (3 corresponding to up or down 
spins, or the form cLcjj, have to be carried out only 
for the active sites. Therefore, the construction of auxil- 
iary operators containing various combinations of one, 
two, three and four fermion operators can be obtained in 
a much faster way. This is especially important for the 
quadratic auxiliary operators since their multiplication 
during the diagonalization of the superblock Hamiltonian 
is the most time consuming part of the MS-DMRG. It 
is worth to mention, that the auxiliary operators gener- 
ated by the left block states have be to calculated for the 
empty sites as well, because during the backward sweep 
the structure of the right block might change and non- 
zero matrix elements in the superblock Hamilton opera- 
tor could be generated. 

This procedure has the advantage, that when doubling 
the number of sites (doing a calculation on a larger ba- 
sis set), the computational time will not increase signifi- 
cantly during the first half sweep and a much larger sub- 
set of block states (larger M™ m ) can be used very effi- 
ciently to generate a better environment block. For opti- 
mal ordering this gives a much faster convergence and due 
to the better environment the crossover between the envi- 
ronment and truncation error can be reached within a few 
steps. After the end of the first half sweep M™" = M r mm 
can be used. 



C. Application of Abelian point-group symmetry 
in QC-DMRG 

When the DMRG procedure is applied in quantum 
chemistry to calculate the energy of molecules the use 
of Abelian point-group symmetry reduces the dimension 
of the Hilbert space to be considered. Before proceed- 
ing with the numerical results we show how to use this 
symmetry in QC-DMRG. The irreducible representation 
of each molecular orbital is generated by standard quan- 
tum chemistry programs, like the MOLPRO [3{| program 
package. The quantum numbers of the irreducible rep- 
resentations are very similar to the momentum quantum 
numbers ki used in MS-DMRG. The Ki quantum num- 
ber operator for an orbital defined on the (0, |, |, IT) basis 
states in QC-DMRG is written as 
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where ki can take values from one to eight according to 
the symmetry of the orbital. The symmetry quantum 
number operator for basis states of two orbitals is deter- 
mined by using the standard character table as 

Kij = T(Ki, Kj), (27) 
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Using Eq. we can determine the symmetry quantum 
numbers of all the left and right block bases states and 
for a given target state we can restrict the number of 
electrons and the total symmetry (Ktot) as 



K TOT =T(K L ,K R ) 



(29) 



where Ktot again can take values between one and 
eight. During the renormalization procedure, the sym- 
metry quantum number operators are renormalized in a 
similar way as the particle number operators and they 
are assign to the new renormalized block states. 



VI. NUMERICAL RESULTS USING 
OPTIMIZED ORDERING AND DEAS 

In order to demonstrate the efficiency of DEAS pro- 
cedure we have studied the 1-D Hubbard chain for var- 
ious fillings, orderings, U values and chain lengths. We 
have found that for small U in the range U ~ 0.1 — 1 
MS-DMRG converges much faster than with the ordi- 
nary initialization procedure. Performing calculations on 
half-filled lattices up to 30 sites the error margin set by 
TRE max ~ 10~ 5 has been reached within one and a half 
or two sweeps. As the strength of U increased the proce- 
dure selected out more block states but the method still 
converged much faster. As an example Fig. ^] shows 
the result obtained for the half-filled, L = 18 chain for 
U = 0.5,1,2,4. In order to obtain an absolute error of 
10~ 4 , TRE max was fixed to 10~ 5 and M min was chosen 
as 400 and we used the ordering shown in the 6th panel 
of Fig. [SJ For [7 = 4 we have cut the maximum num- 
ber of block states at 2500 giving a maximum value of 
the truncation error of the order of 10 -4 . The AS vector 
is AS = [13,5,14,4,12,6,15,3,11,7,16,...] and the HF 
vector is [5, 6, 7, 8, 9, 10, 11, 12, 13]. It can be seen in the 
figure, that in all cases the desired accuracy was obtained 
within two sweeps in contrast to the results obtained by 
Nishimoto et al |6j using the standard initialization pro- 
cedure. It is even more important to emphasize that after 
the first (second) sweep the calculated momentum distri- 
bution (n k ) = (52i t0l c\ a c ia ) for U = 1(U = 4) agreed up 
to 3 digits with the values obtained by real space DMRG 
with periodic boundary condition. 

Similar results were found for the molecules studied. 
Results for the CH2 molecule for fillings (number of elec- 
trons/number of orbitals) 6/13 and 6/23, H2O for 8/24, 
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using DEAS procedure. 



N 2 for 10/26 and F 2 for 18/18 are shown in Fig.[T3 In all 
cases the AS vector was determined in a self-consistent 
way by first running a quick full sweep with some 64 block 
states and then it was constructed based on the decay of 
the site entropy values described in the previous section. 
In all cases we have used the modified Cuthill-McKee or- 
dering, according to the method described in Sec. IV. It 
is clear from the figure that the DMRG with DEAS pro- 
cedure converged very fast and usually 90-99 percent of 
the correlation energy was obtained within the first one 
half sweep. Therefore, the crossover between the environ- 
ment and the truncation error can be reached much faster 
and the accuracy of the method is finally determined by 
TRE max as was shown in Ref. |9j. As an example for 
CH2 for 6/13 filling the same accuracy that was shown 
in Fig. El was reached in 13 iteration steps with DEAS 
procedure using the same parameter sets. For F2 the 
same accuracy that was shown in Fig. 1101 was reached in 
23 iteration steps. 

It is very important to mention, that besides the very 
fast convergence, if an optimal ordering is obtained and 
a good initial condition is used then the number of block 
states can be reduced significantly during the process of 
renormalization and finally a block can describe the to- 
tal system very accurately. As an example we have set 
Mmin = 4 and the energy of the CH2 molecule with 
13 orbitals was obtained after the second sweep up to 
7 digits accuracy with Mi = 21, M r — 4, dim (J = 107 
and for F2 18/18 up to 5 digits accuracy with M/ = 
120, M r = 4,dimC7 = 220 in contrast to the dimension 
of the exact solution which are 230230 and 9075135300 
calculated respectively by the (2L)!/((2£ - N)\Nl) for- 
mula. The same was obtained for the half-filled Hubbard 
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FIG. 15: Percent of correlation energy as a function of iter- 
ation step obtained by DEAS procedure for the CH2, H2O, 
N2, F2 molecules. 



model up to 4 digits accuracy with L — 18, U = 1/2 is 
Mi = 280, M r = 4dim£ = 562. This feature is expected 
to have a close relationship to quantum data compression 
[21II22I . Therefore, it seems to be promising to examine 
the separability of reduced density matrices and their ef- 
fect on the number of block states. A more detailed study 
of the relationship between site entropies and the dimen- 
sion of the Hilbert space of the superblock Hamiltonian 
based on Suchmacher quantum data compression is in 
progress. 



VII. SUMMARY AND FUTURE 
PROSPECTIVES 

We have studied the effect of ordering of lattice sites 
in the momentum space version of the density matrix 
renormalization group (MS-DMRG) method by solving 
the 1-D Hubbard chain and various molecules. This was 
done by calculating site and block entropies, the separa- 
bility and entanglement of the target state. Our findings 
are listed below: 

(1) We have shown that the ki sites or molecular or- 
bitals lying closer to the Fermi surface have larger in- 
formation content and correspond to a larger mixture 
described by larger Neumann entropy and participation 



number. 

(2) We have shown that if all the sites with large en- 
tropy are placed close to one end of the chain then this 
improper ordering can lead to vanishing mutual entropy, 
lack of entanglement and quantum information exchange 
between the system block and environment block. This 
prevents the method to convergence to the target state. 
Therefore, the proposal by Xiang to place the highly cor- 
related sites as close as possible is not a sufficient condi- 
tion alone. 

(3) If lattice sites with large entropy are moved to 
the two ends of the chain the block entropies can be 
maximized giving rise to a larger entanglement, but the 
convergence of the method can still be very slow or the 
method does not converge at all. Numerical results in- 
dicate that if these sites are moved to the center of the 
chain then a very fast convergence can be obtained with 
much smaller subset of block states. We have found that 
the mutual entropy has sharp peaks when the most ac- 
tive sites are just between the system and environment 
blocks when the energy also improves significantly. Fur- 
thermore, in order to have a finite interaction between the 
blocks some sites with smaller but finite entropies have to 
be placed close to the ends of the chain. This keeps the 
block entropy above a critical value for several iteration 
steps. We have shown that these conditions are often not 
satisfied when using the Cuthill-McKee algorithm. 

(4) We have developed a new initialization procedure 
which gives finite block entropies even during the first 
half sweep of MS-DMRG and reduces the volume of sep- 
arable basis states of the target state. This new method 
extends the active space in a dynamical fashion result- 
ing in a very accurate set of renormalized block states 
and a very fast convergence. For the models that have 
been studied 90-99 percent of the correlation energy was 
obtained within the first half sweep. 

(5) We have shown that MS-DMRG is very sensitive 
to the ordering and the initial conditions. A good start- 
ing configuration and an optimized ordering can result in 
very accurate blocks with very limited number of block 
states which still describes the system with the accuracy 
determined in advance. 

We expect that DMRG can be a very good candidate 
for a new method on quantum data compression and on 
quantum error correction |37l l38l l39l | in which case the 
fidelity can be defined in advance by TRE max . This re- 
search is in progress. 

Since MS-DMRG describes a composite system with 
long-range interactions it is expected that the method 
could be used in the context of non-extensive thermo- 
dynamics to study models for which the subadditivity 
of the subsystem entropy does not hold. This would al- 
low one to investigate Tsallis entropy |40j , non-extensive 
mutual entropy and quantum entanglement by using the 
generalized Kullback-Leibler entropy [411 ]. 
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